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Imperfections change essentially the electronic transport properties of graphene. Motivated by a 
recent experiment reporting on the possible application of graphene as junctions, we study transport 
properties in graphene-based junctions with single localized defect. We solve the Dirac-Bogoliubov- 
de-Gennes equation with a single localized defect superconductor-normal(graphene)-superconductor 
(SNS) junction. We consider the properties of tunneling conductance and Josephson current through 
an undoped strip of graphene with heavily doped s-wave superconducting electrodes in the limit 
Idef ^ i ^ We find that spectrum of Andreev bound states are modified in the presence of 
single localized defect in the bulk and the minimum tunneling conductance remains the same. The 
Josephson junction exhibits sign oscillations. 

PACS numbers: 74.45. +c, 74.50. +r, 73.23.Ad, 74.78. Na 



I. INTRODUCTION 



Recent exciting developments in transport experiments 
on graphene has stimulated theoretical studies of super- 
conductivity phenomena in this material, which has been 
recently fabricated [J Q • A number of unusual features 
3 of superconducting state have been predictedjl, [1] 
which are closely related to the Dirac-like spectrum of 
normal state excitation[^-Q. In particular, the uncon- 
ventional normal electron dispertion[Tol| has been shown 
to result in a nontrivial modification of Andreev reflection 
and Andreev bound states in Josephson junctions [llj 
with supercomducting graphene electrodes [l^ [Tsj. 




FIG. 1: A graphene undoped ribbon is contacted by two su- 
perconducting leads. The charge carriers tunnel from one 
lead to another via multiple tunneling states formed in the 
graphene strip. A defect placed inside the strip. 



Other interesting consequences of the existence of 
Dirac-like quasiparticles can be understood by study- 
ing superconductivity [l3| in graphene [T5l - [l9| . It has 
been suggested that superconductivity can be induced 
in graphene layer in the presence of a superconducting 
electrode near it via proximity effect |20l - |22j . 



In this work, we study Josephson effect and find bound 
state in graphenejl^ for tunneling SNS junction with the 
presence of a single localized defect [2J|. In this study, we 
shall concentrate on SNS junction with normal region 
thickness L ^ ^, where ^ is the superconducting coher- 
ence length, and width W which has an applied gate volt- 
age U across the normal region pBi [26j . In the frame of 
the limit Idej <C i ^ ^ cosidered by Kulik we investigate 
tunneling conductance in SNS junctions with presence a 
single localized defect and find that Andreev levels are 
modified, the minimum tunneling conductance remains 
the samefHIH. 



II. JOSEPHSON EFFECT IN SUPERCONDUC- 
TOR/NORMAL (GRAPHENE) /SUPERCONDUCTOR 
JUNCTIONS WITH A SINGLE LOCALIZED 
DEFECT 



We consider a SNS junction with a single localized de- 
fect which is involved in a graphene sheet of width W ly- 
ing in the x — y plane extends from x = —L/2 to x = L/2 
while the superconducting region occupies \x\ > L/2 (see 
Fig. 1). The SNS junctions can then be described by the 
Dirac-Bogoliubov-de-Gennes (DBdG) equations [30|. 



Hs- Ef + U a 

A* Ef~U - Hs 



ips = eips 



Here , ips = {■ipAs,tpBs,^*As^-'^*Bs)^ ^ = -"2, wi, ^2) 
are the 4 component wave functions for the electron and 
hole spinors, the index s denote K 01 K for electrons 
or holes near K and K points, s takes values K{K ) 
for s = K(K ), Ef denotes the Fermi energi, A and B 
denote the two inequivalent sites in the hexagonal lattice 



of graphene, and the Hamiltonian Hs is given by 



-ihvp[ijxdx + sgn{s)ijydy 



(1) 



In Eq. [11 Ui? denotes the Fermi velocity of the quasi- 
particles in graphene and sgn(s) takes values ± for s = 
K{K ). The 2x2 Pauli matrices ai act on the sublattice 
index. The excitation energy e > is measured relative 
to the Fermi level(set at zero). The electrostatic poten- 
tial U and pair patential A have step function profiles, as 
in the case of a semiconductor two-dimensional electron 
gas UMI: A(r) -j. Aq expT **/2 -[/ for x ±oo. We 
assume non-interacting electrons in the normal region, 
therefore, A(r) ee 0, [/ = for |x| < L/2. The reduction 
of the order parameter A(a;) in the superconducting re- 
gion on approaching the SN interface is neglected; i.e., we 
approximate parameter A(x) as we have done it above. 
As discussed by Likharev |34| , this approximation is jus- 
tified if the weak link has length and width much smaller 
than ^. There is no lattice mismatch at the NS interface, 
so the honeycomb lattice of graphene is unperturbed at 
the boundary, the interface is smooth and impurity free. 
Zero magnetic field is assumed. 

Solving the DBdG equations, we gain the wave- 
functions in the superconducting and the normal regions. 
In region I {I II), for the DBdG quasiparticles moving 
along the ±x direction with a transverse momentum 
ky ^ q and energy e, the wave-functions are given by 



^ = exp{iqy + iksX + Kmx) 



/ exp(— im/3) 
exp(i7 — im/3) 
exp(— irn0/2) 
\ exp(j7 — im(f>/2) 



exp{iqy — ik^x -\- Kmx) 



/ exp(i77i/3) 
exp(— 17 -I- im/S) 
exp(— ?m0/2) 
\ exp(— Z7 — im(j)/2) 



The parameters /3, 7, kg, k are defined by 
(3 ~ arccos(e/Ao), 7 — arcsm[hvpq/{Uo + Ep)], 
ks = ^iUo + EF)yihvp)^-q\ K = (t/o + 
i?i?)Ao sin(/3)/(?i2wpA:s) and m ~ ± denotes region 
I {1 1 1), m = + for / and m = — for /// correspondently. 
Further we assumed that the Fermi wave length in 
the superconducting region much smaller than the wave 
length Af in the normal region and Uq ^ Ep,e. Since 
kl < Ep/hvp, this regime of a heavily doped supercon- 
ductor corresponds to the limits 7 — >■ 0, kg Uq/Hvp, 
K (Ao/hvp) sin(;3). 

Now we analyze the spectral properties of a graphene 
ring. Formally we find solutions for graphene ring and 
then extending ring in the scale, match external boarder 
of ring with superconducting regions of junction and 
fix internal part implying which as defect. For that we 
devide the region // into three areas: a, 6, c (see Fig. 
2). We solve the DBdG equations for the area 5, while 
area c is the area where placed defect and area a would 




FIG. 2: The normal region II is devided by three areas a, b, c. 
The DBdG equations are solved for area b and determined by 
the "infinite mass" boundary conditions induced by V{r) — >■ 
+00 in a and c areas correspondently. 



be extended and matched with superconducting regions. 
The two valleys s± decouple, and we can solve equa- 
tions separately for each valley, Hgijjs = (e + sEp)'ips, 
Hs = Hq + sV{r)az- The term proportional to Uz in 
Hamiltonian is a mass term confining the Dirac electrons 
in the area b. Rewrite the Hamiltonian in the cylindrical 
coordinates and since Hg commutes with = Iz + ■^o'z, 
its electron-eigenspinors 4'e are eigenstates of Jz [1^, 



*e(r,a) 



exp{id{n ~ l/2)a) J<i(„_i/2) (fc(e)r) 
exp{id{n + l/2)a) Jd(n+i/2)(fc(e)'') 



with eigenvalues n, where n is a half-odd integer, n — 
di, d|, . . . and Jd{n-i/2){k{€)r) is the Bessel function of 
(n — 1/2) order. In the x ~ y plane d denotes moving 
direction of correspondent quasiparticle, d = + for the 
quasiparticle moving toward x = L/2 and d = — for the 
quasiparticle moving toward x = —L/2 direction corre- 
spondently. Further we are interested in to find zero en- 
ergy states |22|. In this case the DBdG equations posses a 
general symmetry with respect to the change in the sign 
of energy [sll, 



(2) 



where we denote u = (ui, U2) and v ~ (wi, V2). Thus, for 
a set of zero modes {ui,Vi) enumerated by a certain index 
i we should have. 



V. 



(3) 



In the same manner as for electron hole-spinors have 
view, 



*/i(r-,a) = 



-exp(zd(n-|- l/2)a')Jd(«+i/2)(fc'(e)?') 
exp{id{n - l/2)a' )Jd(n^i/2) {k' (e)r) 
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with the definitions 



a(e) = a,rcsm[hvFq/ + Ep)], 
a (e) — arcsin[/iWi?(7/(e — Ep)], 

fc(e) — {hvF)^^{e + Ep) cos(a), 
k (e) = (?iw_F)^"^(e — Ep) cos(a), 



(4) 
(5) 

(6) 
(7) 



The angle a £ (— 7r/2,7r/2) is the angle of incidence 
of the electron (having longitudinal wave vector fc), and 
a is the reflection angle of the hole (having longitudinal 
wave vector k ) [stI . [38| . To obtain an analytical approx- 
imation of the spectrum, we use the asymptotic form of 
the Bessel functions for large r. This indeed is the de- 
sired limit as rfc(e) « rdefk{e) oc rdef/L <^ 1, where rdef 
is the defect radius (the radius of the area c) and deter- 
mine for all eigenvalues n = di. In this limit we impose 
the "infinite mass" boundary conditions at y = 0, W^, for 
which g„ = {n+l/2)TT /W in the area b with V{r) +00 
in the a and c areas consequently (see Fig. 3). Half-odd 
integer values n reflect the tt Berry's phase of closed size 
of a single localized defect in graphene. 

To obtain the subgap (e < Aq) Andreev bound states, 
we now impose the boundary conditions at the graphene. 
The wave-functions in the superconducting and normal 
regions can be constructed as 



(8) 
(9) 



where ai (61), 02(62) are the amplitudes of right and left 
moving DBdG quasiparticles in region I(III) and a(6) 
and c{d) are the amplitudes of right(left) moving elec- 
trons and holes, respectively, in the normal region[l^. 
These wave functions must satisfy the boundary condi- 
tions, 

^l\x=-L/2 = *//U=-L/2, «'//U=L/2 = «'///U=L/2-(10) 

Since the wave vector ky parallel to the NS interface and 
different wave vectors in the y-direction are not coupled, 
we may solve the problem for a given ky = q and we 
can consider each transverse mode separately. To leading 
order in the small parameter AqL/Hvp we may substitute 
a{a ) a(0), k{e){k (e)) — )> fc(0). After some algebra we 
obtain equation 



cos(2/3) 
-sin(2/3) 



sin^(fci) - cos^ {kL) sin^(a) 



cos2(fcL) cos2(a) 
cos(fcL) sin(fcL) sin(a 



1 



C0S((/)) 



(11) 



cos^ (kL) cos2(a) — 1 

Eq. [TT] differs from equation obtained for SNS junction 
without a single defect. Eliminating second term in Eq. 
[TT]we could immediately yield reduction of the equation 
and it earns a essential form for weak SNS junctions [23|. 
The solution of Eq. [TT]is a single bound state per mode, 



A01 



■(-2CA2 -B^+ BV1-4CA2(C + 1)) 




FIG. 3: Tunneling conductance (normalized per mode) of the 
graphene SNS junctions with a single localized defect versus 
Fermi energy, calculated from Eq. 1141 The tunneling conduc- 
tance exhibits oscillatory behaviour. 



where 



A = 



sin^(fc„L) — cos2(fc„I/) sin^(Q;) 
cos2(fc„L) cos2(a) — 1 
cos(A:„L) sin(fc„L) sin(a) 
cos2(fc„L) cos2(a) — 1 



B 



D 



C=i + l(i-sin^(^)), 

cos^{knL) cos^{a) — 1 
cos2(fc„L) cos2(q!) — cos(2/c„L) ' 



(12) 
(13) 

(14) 
(15) 



We don't have a simple analytic expression for the 
0-dependance but we obtained modified Andreev lev- 
els with the presence a single localized defect in the 
bulk. The conductance of the graphene strip is expressed 
through the transmission probability by the Landauer 
formula. 



go = 4eV/i, 



(16) 



where n{fj.) >• 1 is given by n(^)=Int(fc„M^/7r + 1/2). 
Substitution transmission probability into Eq. [16] gives 
the conductance (normalized per mode) versus Fermi en- 
ergy (see in Fig. 3). 



III. JOSEPHSON CURRENT 

The Josephson current at zero temperature is given by 

^('^) = -Y^y^ deY,Pnie,<ly)e (17) 



r!=0 



where the factor of 4 accounts for the twofold spin and 
valley degeneracies. Substitution of pn{s,(l)) = ^[^ ~ 
e„(0)] into Eq. [T7] gives the supercurrent due to the 
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FIG. 4: Josephson current Ic of graphene-based SNS junc- 
tion with single localized defect (length L short compared to 
the width W and superconducting coherence length as a 
function of the Fermi energy ^ in the normal region. 



discrete spectrum 



1 



An sin ( 



Contributions to the supercurrcnt from the continuous 
spectrum are smaller by a factor L/^ and may be ne- 
glected in the short-junction regime [6]. For L -C 
the summation over n may be replaced by an integra- 



tion. The resulting Josephson current Ic upon substitu- 
tion (/) — > 7r/2 is plotted as a function of fi in Fig. U) 



IV. CONCLUSION 

In conclusion, we have shown that a Josephson junc- 
tion in graphene can carry a nonzero supercurrent even if 
the Fermi level is tuned to the point of zero carrier con- 
centration. At this Dirac point, the current-phase rela- 
tionship has a resonant behaviour due to single localized 
defect. This unusual " quasidiffusive" behaviour of the 
Josephson effect in undoped graphene should be observ- 
able in submicrometer scale junctions. It is found that 
the current has a peak-like structure, with alternating 
signs of the peaks. The resulting nonequilibrium Joseph- 
son current for a given phase difference thus oscillates as 
a function of chemical potential, i.e the junction displays 
so-called 7r-behaviour. 
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